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ABSTRACT 

We discuss the possibility that vibrationally 
excited H ^ may be collisionally de-excited, so providing 
a heating mechanism for interstellar clouds which 
operates by coupling the stellar radiation to the gas. 

The majority of excitations in the Lyman and Werner 
bands of Hg return the molecules to the ground electronic 
state in a vibrationally excited level, the most 
favored being v" = 7. We compare the heating rate 
obtained in this way with other mechanisms which have 
been postulated, and present the results of calculations 
of temperature as a function of depth into clouds of 
different densities. It appears that this mechanism 
is a significant one, which should be taken into account 
in detailed models of dense clouds. 



I . INTRODUCTION 


There has been a continuing effort to understand the 
heating and cooling processes in the interstellar medium 
since Spitzer’s classic papers (1948-50). The most recent 
review is that of Dalgarno and McCray (1972) who list a very 
large number of collisional processes contributing to cooling 
and discuss time-independent and time-dependent heating 
effects. This article presents a heating mechanism not 
considered hitherto, which may be applicable to dense 
interstellar Hg clouds. It can effectively couple stellar 
radiation ~10008 with the gas by the collisional de-excitation 
of vibrationally excited Hg. 


II. THE MECHANISM 


It has been shown by Stecher and Williams (1967) that 
a significant fraction of all excitations in the Hg (v M = 0) 

Lyman bands by the absorption of ~100()8 radiation results in 
destruction of the molecules because the cascade is into the 
vibrational continuum of the ground state. Accurate calculations 
by Dalgarno and Stephens (1970) show this fraction to be 23%. 

The remaining 77% of the Lyman band excitations, and all the 
Werner band excitations (for \>9128) return the molecules to 
vibrationally excited bound levels in the electronic ground 
state. The initial distribution among these levels is known 
(Allison and Dalgarno 1969, Cartwright and Drapatz 1970), v" 
in the range 4-10 being most favored. Since Hg has no dipole 
moment, transitions between the vibrational levels of the 
ground state are forbidden. The vibrational transition 

probabilities are therefore very small; for example, the 

-7 -1 

v" = 1 -* v" = 0 transition has a rate 8.56 x 10 sec 
(Dalgarno, private, communication, 1979 . These values imply 
long radiative lifetimes, and so it is of interest to assess 
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the importance of collisional de-excitation of these levels. 

In such collisions, the vibrational energy removed is trans- 
formed into kinetic energy. 

Formally, we may treat the problem by writing down the 
equation for the number density of E^ in any particular 
vibrational level in terms of absorption of radiation, and 
of relaxation by radiative and collisional processes. In 
equilibrium, there are fifteen coupled algebraic equations, from 
which the mean level populations may be determined. Then the 
heating function is given by the relaxation rates times the 
appropriate energy transfer. 

Thanks to the work of Dalgarno and his associates much 
of the necessary data concerned with radiative processes is 
available. However, the collisional data is extremely poorly 
known, if at all. Most of that which exists refers only to 
the 0-1 band. Consequently, with the present level of knowledge 
the full set of equations cannot be solved. However, we assert 
that useful information can be obtained from a two- level system 
in which the lower level represents v" = 0, and the upper 
level represents all 0 < v" ^ 14. We shall use a radiative rate 
for the upper level corresponding to the inverse of the radiative 
lifetime of the v" = 7 level and a mean energy released per 
relaxation collision equivalent to this level also. This assumes 
that relaxation is complete when it occurs. 

III. THE TWO- LEVEL SYSTEM 

For this system instead of fifteen equations, we obtain 

^ nCH^) = agn (H 2 ) - a* p* n (H 2 > - An(H*) - Rn(H*) (1) 

where nCH^) represents the number density of vibrationally 
excited Hg molecules and nCHg) that of H v" = 0 molecules. 

The rate of excitation from the X, v" = 0 state is P, and 

the fraction of all excitations returning to an X, v" > 0 bound 
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state is a. In a similar way, a*p* represents the fraction 
of all excitations of H 2 not returning as H 2 « This term- 
includes direct photodissociation and photoionization 

>{c 

H 2 which can occur in the interstellar medium (Stecher and 
Williams 1967) . The symbols A and R represent the radiative 
relaxation (Einstein A-value) and the collisional relaxation 
(H “* K^) processes, respectively. Equation (1) assumes that 
nCHg) is constant in time, i.e. that the production and 
destruction rates for H 2 are equal. The mechanism may also 
be valid in non-steady situations. In equilibrium 

n(H*)/n(H 2 ) - 00/ (A + R + a*0*) (2) 

and the heating rate is now 

e o n(H*) R 

i.e. 

= n( - H 2^ e o {a + R + 

where e Q is the mean energy released in the relaxation. It 
is convenient to discuss the term in the curly brackets 
separately from the rest of equation (3). Therefore, we 
introduce the notation 

§ = R/(A + R + a*p*) (4) 

for the factor containing all the relaxation terms. This 
represents the coupling between the radiation energy and 
heat energy in the gas. 

IV * RELAXATION PROCESSES 

(a) Electrons - Vibrational excitation in slow e-Hg collisions 
has been studied experimentally (Ehrhardt et al. 1968, Trajmar 
et al. 1969, and in particular, near the v" = 1 vibrational 
excitation threshold, Burrow and Schulz 1969) and also 
theoretically (Henry, 1970). From the relations between 
the excitation and de-excitation cross-sections 
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(Seaton 1958) we find the v" = 1 -* v" = 0 de-excitation cross- 

- 17 2 

section is ~2 x 10 cm for any thermal energy which is much 
less than the vibrational energy. We therefore, adopt for 
electrons : 

Xa u,,S _, n -10 +3-l /cN 

coll coly i -♦ o e ~ cm s , (5) 

while recognizing that the corresponding rate for higher 
v" -» v" - 1 transitions may be faster. 

(b) Atoms and Molecules - Takayanagi (1963) has reviewed the 

subject of collisionally induced vibrational transitions 

and presents criteria which make it clear that - under conditions 

appropriate to the present discussions - such transitions are 

unlikely. The existing evidence supports this view. 

Relaxation of H 0 in collision with H 0 is very slow; measured 

Z -j^g z 

times are on the order of 10 /n^g) seconds at room temperature 

and even longer at lower temperatures (Kiefer and Lutz, 1966; 
DeMartini and Ducuing, 1966). Relaxation of E^ by atoms H 
and He has been examined theoretically. For example, Micha 
and Rotenberg (1970) have considered He, Hg collisions. 

They find cross-sections which are negligible for our purpose. 

The case of H-atom relaxation is made slightly more favorable 

jjg 

by the possibility of exchange, H + E^ H 2 + H., enabled by 

tunnelling through the appropriate barrier. Truhlar and Kupperman 

(1970) have used an appropriate Hg surface and find that the 

probability of a non-adiabatic collision is greater than 1/2 

for a relative translational energy of some tenths of an eV. 

Allowing a reasonable barrier shape, we find that a cross- 
-19 2 

section of ~10 cm would be appropriate under the relevant 
conditions. We conclude that such interactions are unimportant 
compared to electrons. But it is possible that in some 
extreme situations they may be important. We, therefore, present 

the number densities of some de-exciting agents at which 

-7 -1 

collisional and radiative de-excitation rates (~10 s ) are 

10 7 

equal. For H^, H, e, p these number densities are ~10 , 10 , 
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10 and 10 cm , respectively. 

V. DISCUSSION OF THE FACTOR 5 
From the vibrational transition probabilities supplied by 

A. Dalgarno (private communication 1972) , the mean cascade 

-7 -1 

rate A which we should insert in (4) is 3 x 10 s . From 
(5) we adopt for R 

R = 10“ 10 n(e) (6) 

The other term in § is a (3 . The information required to 

sj c 

calculate a , the fraction of excitations of H 0 not returning 
to , is available (Allison and Dalgarno 1969, Dalgarno and 
Stephens 1970, Cartwright and Drapatz 1970) and we present in 
Table I the fractions for each vibrational level, for the 
Lyman and V/erner bands, including only those transitions for 
which \>912S. We choose a mean a* to be 0.2. The rate of 
excitation in the Lyman and Werner bands from H 2> f3 , will be 
larger than (3 since (with due allowance for rotational branching) 
42 times as many levels are involved. The amount of f-value 
associated with each v" level is approximately the same. We 
shall assume f3 z 42j3. In the optically thin case, the lifetime 

of an H 0 molecule in a typical interstellar radiation field is 

10 " _ 
~10 s (Stecher and Williams 1967). This means that a|3 2 4 x 10 

— 1 * ale 

s and from this we calculate that a |3* in the optically thin 

-9 -1 

case is ~ 4 x 10 s . Allowance for direct photodissociation 

* -9-1 

of H 0 increases this to ~ 5 x 10 s . In most situations, 

^ _4 \ 

therefore, we may write § = 3 x 10 n(e) and the difficulty 
lies in the choice of n(e) for the situation considered (Pottasch 
1972) . 


VI. THE HEATING RATE 

The heating rate is now given by 

= n(H 2 ) ape Q § (7) 

in which extinction and line absorption problems 
are all contained in (3. We calculate extinction due to the 
grains at a rate of 3 magnitudes per kiloparsec of the mean 
interstellar medium (Stecher 1SS9) , assuming a constant gas/ 

- 5 - 



dust ratio. There is also competition for the ~ 1000$ radiation 

by the photoionization of carbon, which has a threshold at 

- 17 2 

11.1 eV and a nearly constant cross section of 1.1 x 10 cm 

over the first few eV of absorption (Bates and Seaton 1949). 

-Werner (1970) has shown how the fractional abundance of carbon 

varies with depth inside clouds of various densities on the 

assumption that the electron density arises solely from the 

photoionization of carbon. This assumption may be in error. 

The line widths as a function of column density are from 

-12 

Cartwright and Drapatz (1970). We choose = 5 x 10 erg. 

In any given situation we may now evaluate#^. 

VII. OTHER HEATING MECHANISMS 

Besides the familiar suggestions of cosmic rays and X-rays 

(see e.g. Dalgarno and McCray 1972) as heat sources for the 

interstellar medium, and the recent idea of Watson (1972) 

concerning grains, there are a number of other heating mechanisms 

which may be significant. These have in common with the Hg heating 

rate the ability to couple stellar radiation energy to the 

gas as distinct from cosmic- or X-ray energy. They depend on 

photodissociation and photoionization of molecules and ions. 

In general, each such event will release several eV into the 

gas. This process may be so efficient at using relatively 

low energy continuum radiation, that the rate-limiting 

mechanism is usually the formation rate of the relevant species. 

The photodissociation of a diatomic molecule XY of 

dissociation energy D by radiation of energy E gives a heat 

gain (E-D) . The fastest production rates via radiative 

association are ~ 10" cm s~ (Julienne and Krauss 1972) and 

thus the greatest heating rate that can be obtained is 
—28 —3 —1 

~ 10 n(X)n(Y)erg cm s. Grain surface reactions may improve 
this rate somewhat. This is unlikely to be significant in 
dense clouds. 

Photodetachment in negative ions also provides a possible 
heat source. Using the radiative attachment rate of de Jong (1972) 


6 



for H+e we find a heating rate ~ 10 ®®T n(H)n(e)erg cra - ^s - ^ which 

clearly has little significance in dense clouds but may be 

important in the intercloud region. 

In section III we emphasized that allowance should be made 

for- direct photodissociation of E ^ (v s 3). Each such 

photodissociation supplies heat and using III we estimate a 

-24 — 3 - 1 

heating rate of ~ 10 nCE^) erg cm s in the optically 
thin case at the edge of clouds. In the interior this rate 
falls off quickly since two photons are required to make it 
work. But this rate may be significant in low density clouds 
where a high H ^ production rate exists. 

VIII. DISCUSSION OF THE SIGNIFICANCE OF THE HEATING RATE 

To determine the significance of we compare it 

with other heating rates which have been suggested. All of 

these heating rates fall off with depth. There is considerable 

\ 

doubt that the lower energy cosmic rays penetrate denser 

clouds at all (Solomon and Werner 1971). X-rays of 1/4 keV 

are absorbed with an optical depth of unity in a column density 
19-2 

of 6 x 10 cm . The grain mechanism of Watson competes 
directly for the photons operating the X^CEL) mechanism. Because 
the modes of decrease are complicated, we first show in Table II 
heating rates at zero depth into clouds for the various processes. 

The rate «/r*( H„) is given for an incident flux of 10 photons cm" s" 
the optically thin case, and on the assumption that n(e) = n(C + ), 
a minumum case. Clearly, these zero depth rates are not 
very meaningful in an absolute sense but this table shows 
that-especially in denser regions - the heating rate ^(HL) 
must be considered. 

A clearer and perhaps more realistic way to assess the 
importance of ^(Hg) is to couple it with a cooling rate to 
obtain a temperature at various depths in a cloud ignoring 
all other heating mechanisms. We present the results of 
these calculations in Table III. The incident radiation from 
at ~ 1000& was assumed to be 10® photons cm“^s“ ^ . 
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The cooling rates are taken from Dalgarno and McCray (1972) 
and include cooling due to excitation in C, C + , O, CO, CN 
and H 2 by e and by Hg. However, the rates quoted must be 
modified for the high density regions considered here due to 
collisional de-excitation of these species [cf. Bahcall and 
Wolf 1968]. In the calculation of vP(H 2 ) the equivalent 
widths of Cartwright and Drapatz (1970) were used, and allowance 
was made for extinction and the photoioniza# ioh of carbon, the 
electron density being assumed to arise solely from this process. 

Thus we have a minimum Jv(H 2 ) ignoring possible electrons from 
other sources. We also include temperatures which would be 
obtained if n(e) were 30 times greater. Finally, we show 
the temperatures which would be obtained if heavy atoms were 
depleted by sticking to grains in the manner suggested by 
Field, Goldsmith, and Habing (1969) . The depletion factor is 
10 times, and the effect is to reduce the cooling rate. The 
temperatures are given at various depths within the 
cloud measured from the edge, the incident flux of radiation 
being chosen as 10^ photons cm - ^s”‘*'8 - ^ in the relevant wave- 
length range. Note that these temperatures are not precise 
calculations, but are a measure of the importance of jr£(H 2 ). 

It is clear from the cases considered that the heating 
from H 2 is confined to the edges of clouds, as is the heating 
from any other suggested mechanisms. The central temperatures 
of dense regions must therefore be very low indeed. The 
relatively high temperatures at the edges of clouds will 
provide a confining pressure for clouds which are not gravitationally 
bound. The temperatures obtained are not inconsistent with 
observations (Heiles 1971) and this fact supports the view that 
heating by H 2 is an important mechanism in dense clouds. 
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TABLE I 


>}c 

a „ represents the fraction of all excitations 

in the Lyman and V/erncr Lands not returning to 

>[< 

groundstate E^ , for each level v" 


v u 

* 

a v „ (Lyman) 

a v „(7erner) 

v" 

* 

a v „ (Lyman) 

>jc 

a yM (Werner) 

1 

0.19 

0 

8 

0.16 

0.05 

2 

0.21 

0.03 

9 

0.17 

0.03 

3 

0.16 

0.01 

10 ' 

0.17 

0.02 

4 

0.17 

0.02 

11 

0.15 

0.05 

5 

0.18 

0.06 

12 

0.14 

0.08 

6 

0.18 

0.02 

13 

0.18 

0.15 

7 

0.16 

0.02 

14 

0.34 

0.51 


a v „ is calculated from data in Allison and Dalgarno (1969), 
Dalgarno and Stephens (1970), and Cartwright and Drapatz (1970). 
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TABLE II 


Approximate Heating Rates at Zero 

-3 - 1 

Depth into Eg Clouds (erg cm s ) 


<h£) 


2 x 10“ 27 n 2 (Kg) 


Cosmic Rays 


5 x 10“ 27 n(H + Hg) 


X-Rays 


2 x 10“ 26 n(H + Hg) 


Photodetachment from grains 


10“ 26 n(H + Hg) 


-33 

Photodestruction of XH, maximum 10 n(H + Hg) 
(formed in radiative association) “ 


s-35 2 


Photodestruction of XY, maximum 10“ n (H + Hg) 
(formed in radiative association) ^ 


Photodetachment of H“ 


10“ 30 T n(H)n(e) 


Photodissociation of Hg(v" s 3) 


10" 24 n(Hg) 


(but see discussion in VII) 



TABLE III 


Temperatures at Selected Depth Within H 0 Clouds of 
Different Densities Heated Solely by H 2 

2 

Depth(cm) n(H 2 ) = 10 Depletion No Depletion 


1 = io 18 

high n(e) 

19 

14 


low n(e) 

12 

8 

1 = 10 19 

high n(e) 

15 

11 


low n(e) 

9 

«10 


n(H 2 ) = 10 3 



1 = io 17 

high n(e) 

25 

14 


low n(e) 

13 

8 

1 = io 18 

high n(e) 

12 

8 


low n(e) 

«10 

«10 


n(H 2 ) = 10 4 



1 = 10 16 

high n(e) 

65 

24 


low n(e) 

18 

13 

1 = io 17 

high n(e) 

«10 

«10 


low n(e) 

«10 

«10 


Temperatures calculated to be more than a few degrees less than 
10°K are written «10°K, since the cooling curve is not 
particularly accurate- at these temperatures. 
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